Precise Predictions for Z + 4 Jets at Hadron Colliders 
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We present the cross section for production of a Z boson in association with four jets at the Large 
Hadron Collider, at next-to-leading order in the QCD coupling. When the Z decays to neutrinos, 
this process is a key irreducible background to many searches for new physics. Its computation 
has been made feasible through the development of the on-shell approach to perturbative quantum 
field theory. We present the total cross section for pp collisions at y/s = 7 TeV, after folding in the 
decay of the Z boson, or virtual photon, to a charged-lepton pair. We also provide distributions 
of the transverse momenta of the four jets, and we compare cross sections and distributions to the 
corresponding ones for the production of a W boson with accompanying jets. 

PACS numbers: 12.38.-t, 12.38.Bx, 13.87.-a, 14.70. Hp 



The Large Hadron Collider (LHC) is currently ex- 
tending the energy frontier into uncharted territory, in 
the quest to identify new physics beyond the Standard 
Model of particle physics. Many signals of new physics, 
especially those containing dark matter candidates, lie 
in broad distributions with significant Standard Model 
backgrounds. A first-principles understanding of these 
backgrounds is provided by quantum chromodynamics 
(QCD) and the QCD-improved parton model. The lead- 
ing perturbative order (LO) in the QCD coupling a s 
gives a good qualitative prediction. Quantitatively re- 
liable predictions require, at the least, next-to-leading- 
order (NLO) accuracy in the QCD coupling. For pro- 
cesses at a hadron collider with many-jet final states, 
NLO computations have long been a formidable challenge 
to particle theorists. 

In this article we present the first NLO QCD results 
for Z boson production in association with four jets at a 
hadron collider, specifically at the LHC. We fold in the 
decay of the Z boson to an e + e~ pair (or equivalcntly 
and include contributions from virtual-photon 
exchange (collectively denoted by Z,j*). This process, 
containing identifiable charged leptons, is a benchmark 
for the closely related process in which the Z decays into 
neutrinos, which appear as missing transverse energy. 
The Z — ¥ vv decay mode generates a key background 
process in the search for supersymmetry, as well as for 
other models that lead to dark-matter particle produc- 
tion at the end of a cascade of strongly-produced new 
particles. Fig. Q] shows a typical signal process, leading 
to the same signature of missing transverse energy with 
four jets and no sharp resonance. We note that another 
approach to estimating this process — combining a mea- 
surement of prompt-photon production with a theoretical 
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FIG. 1: Gluino pair production illustrates a typical signature 
of new physics scenarios: four jets plus a pair of lightest super- 
symmetric particles (LSPs) that escape the detector, yielding 
missing transverse energy. 
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FIG. 2: Sample diagrams for the seven-point loop amplitudes 
for qg — > Zqggg and qQ — > ZqQ'Q'Q, followed by Z — > e + e~ . 
There are also small contributions where the Z boson is re- 
placed by a photon. This process is very similar theoretically 
to the case Z — > vv with missing transverse energy. 



estimate of the .Z-to-photon ratio [Il-Q — also benefits 
from NLO cross sections 

Recent years have witnessed a growing number of 
NLO QCD results using both traditional and on-shell 
approaches [B|- 11 . On-shell methods 12 - lij exploit the 
analytic properties that all scattering amplitudes must 
satisfy, and generate new amplitudes from previously- 
computed ones. Computationally, they scale modestly 
with increasing numbers of external partons. We used 
these methods to compute the production of a W or Z 
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no. jets 


Z LO 


Z NLO 


Z/W + LO 


Z/W + NLO 


Znl{n-\) LO 


Zn/(n-l) NLO 





323.1(0.1)1^1 


428.6(0.3)14 1 


0.1209(0.0001) 


0.1306(0.0003) 
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66.69(0.04)iy9 


82.1(0.1)1^ 


0.1674(0.0002) 


0.166(0.001) 


0.2064(0.0001) 


0.1915(0.0004) 
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19.10(0.02)1^ 


20.25(0.07)l?;o2 


0.1636(0.0003) 


0.166(0.002) 


0.2864(0.0003) 


0.247(0.001) 
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4.76(0.01)1?;^ 


4.73(0.03)1°;^ 


0.1634(0.0004) 


0.169(0.002) 


0.2494(0.0004) 


0.234(0.002) 
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1.116(0.002)1°;** 


1.06(0.01)1°;?$ 


0.1618(0.0003) 


0.172(0.002) 


0.2343(0.0005) 


0.223(0.002) 



TABLE I: Total cross sections in pb for Z, 7* + n-jet production at the LHC, using the anti-fey jet algorithm with R = 0.5. The 
NLO result for Z, 7* + 4 jets uses the leading-color virtual approximation. The fourth and fifth columns give the cross-section 
ratios for Z,j* — > e + e~ production to W + — > e + v production. The final two columns give the ratios of the cross section for 
the given process to that with one fewer jet. The numerical integration uncertainty is in parentheses. The scale dependence is 
given in superscripts and subscripts. 



boson in association with three jets at NLO 0, [j|. The 
predictions are generally in very good agreement with 
data from the Tevatron [HI, E3- (Earlier NLO results for 
W + 3 jets were based on similar techniques and used 
various leading-color approximations [a Bill.) We have 



also calculated [ll] W + 4-jet production at the LHC, 



making use of a leading-color approximation for the vir- 
tual terms that is known to be valid to about 3% for up 
to three associated jets @, ■ We will use the same ap- 
proximation for Z + 4-jet production. The Z + 4-jet pro- 
duction computation is significantly more complex than 
that for W + 4-jet production, because the quark flavor 
structure leads to more partonic subprocesses, especially 
those containing identical fermion pairs. 

We use the same basic setup as in our earlier work [7j, 
9]. Virtual contributions are evaluated via the Black- 
Hat package [l8[, an implementation of on-shell meth- 
ods. We incorporate a number of significant improve- 
ments in automating the assembly of subprocesses and 
in ensuring the numerical stability of the virtual correc- 
tions. To minimize the amount of higher-precision recom- 
putation at points for which an instability is detected, 
only the unstable part, rather than the whole matrix el- 
ement, is recomputed @, Representative virtual di- 
agrams are shown in fig. [2] We include all subprocesses, 
and make the leading-color virtual approximation only in 
W, Z/^f* + 4-jet production. As in ref. [§}, we drop small 
axial and vector loop contributions, along with the small 
effects of top quarks in the loop. 

The remaining NLO ingredients, the real-emission 
and dipole-subtraction terms [20j ]. are computed using 
AMEGIC++ [HI, which is part of the SHERPA pack- 



age [22|. Here we retain the full color dependence. The 



SHERPA-based phase-space integration exploits QCD 



antenna structures [23J, l24j . BlackHat supplies the 
real-emission tree amplitudes, using on-shell recursion re- 
lations 13j and efficient analytic forms extracted from 
Af = 4 super- Yang-Mills theory [25j. We have validated 
the code extensively. Previously, we compared many re- 
sults against MCFM [H for W, Z + 2-jet production. 
Cross sections and distributions at LO suffer from 



strong sensitivity to the unphysical renormalization scale 
Hn and factorization scale /i p entering a s and the parton 
distributions. This dependence is reduced at NLO. This 
issue is especially important at the LHC because of the 
wide range of kinematics probed. This wide range also 
obliges us to choose an event- by- event scale characteris- 
tic of the kinematics when we compute distributions. We 
choose fin = fiF = M = Hip/ 2 as our central scale [11], 



where H' T 



E z . The sum runs over all final- 



state partons i, and E, z ~ 



; y "z + ipf e ~) 2 ; M z is fixed 
to its on-shell value. We follow standard procedure to 
assess scale dependence, varying the central scale up and 
down by a factor of two to construct scale-dependence 
bands, taking the minimum and maximum of any observ- 
able evaluated at five values: /1 x (1/2, 1/V2, 1, -y/2, 2). 

The fixed-order perturbative expansion may break 
down in special kinematic regions, where large logarithms 
of ratios of physical scales emerge. Threshold logarithms 
can affect production at very large partonic center-of- 
mass energies. However, in ref. [ll| it was argued, using 
results for inclusive single-jet production j27|, that at the 
mass scales probed in W, Z/7* + 4-jet production, such 
logarithms should remain quite modest. Tighter cuts can 
isolate regions subject to potentially large logarithms of 
either QCD or electroweak origin. In particular, cuts that 
force the vector boson to large px , the desired region for 
many searches for supersymmetry or dark-matter parti- 
cles, can induce large electroweak Sudakov logarithms. 

In our study, we consider the inclusive process pp — > 
Z + 4 jets at an LHC center-of-mass energy of y/s = 7 
TeV. We incorporate the full Z, 7* Breit-Wigner reso- 
nance and decay the intermediate boson into an electron- 
positron pair at the amplitude level, retaining all spin 
correlations. We impose the following cuts on the trans- 
verse momenta px, and pseudorapidities ?y: p T > 20 
GeV, |rf I < 2.5, > 25 GeV, |?7 jct | < 3, and 66 
GeV < M e+e - < 116 GeV. The lower cut on the lepton- 
pair invariant mass M e + e - eliminates the large contri- 
bution from the photon pole. Jets are defined using the 
infrared-safe anti-fcT algorithm [28} adopted by the LHC 
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FIG. 3: A comparison of the pr distributions of the leading four jets in Z,j* + 4-jet production at the LHC. In the upper 
panels the NLO distribution is the solid (black) histogram and the LO predictions are shown as dashed (blue) lines. The thin 
vertical line in the center of each bin (where visible) gives its numerical (Monte Carlo) integration error. The middle panels 
show the LO distribution and LO and NLO scale-dependence bands normalized to the central NLO prediction. The bands are 
shaded (gray) for NLO and cross-hatched (brown) for LO. In the bottom panel, the dotted (red) line is the LO Z/W~ ratio, 
the dot-longer-dash (cyan) line the NLO Z/W~ ratio, the dot-shorter-dash (brown) line the LO Z/W ratio and the solid 
(green) line the NLO Z/W + ratio. 



experiments. Here we present results for size parameter 
R = 0.5. We order the jets in pt- In comparisons to W- 
boson cross sections we follow exactly the cuts of ref. 11] ; 
the jet cuts are identical. We use the CTEQ6M par- 
ton distribution functions at NLO, and the CTEQ6L1 
set at LO. Electroweak boson masses and couplings are 
chosen as in refs. @, H|. We also use the SHERPA six- 
flavor implementation of a s (/x) and the value of a s (Mz) 
provided by CTEQ. 

In table U we give LO and NLO parton-level inclusive 
cross sections for e + e~ production via a Z, 7* boson, and 
accompanied by zero through four jets. The NLO results 
exhibit a markedly reduced scale dependence compared 
to LO; the improvement becomes stronger as the num- 
ber of jets increases. We also display the ratios of the 
Z to W + cross sections, and the "jet-production" ratios 
of Z + n-jet to Z + (n— l)-jet cross sections. Ratios to 
VF~-boson cross sections can be obtained using the re- 
sults of ref. [ll| . Both kinds of ratios should be less sen- 
sitive to theoretical systematics than the absolute cross 
sections. Indeed, the Z/W ratios show relatively little 
difference between LO and NLO. This ratio changes very 
little under correlated variations of /1 in numerator and 
denominator; hence we do not exhibit such scale varia- 
tion. Varying the R parameter in the jet algorithm, we 



find very similar behavior as in the W case 11 1. 

It has generally been expected that the jet-production 
ratio is roughly independent of the number of jets [3fj| . 
Other than the Z+ l-jct/Z+ 0-jet ratio, which is smaller 
because of the restricted kinematics of the leading con- 
tribution to Z + 0-jet production, the results shown in 
table U are consistent with this expectation. The ratios 
are, however, rather sensitive to the experimental cuts: 
for example, imposing large vector-boson pt cuts makes 
them depend strongly on the number of jets 

In fig. [31 we show the pr distributions of the leading 
four jets in Z, 7* + 4-jet production at LO and NLO. The 
predictions are normalized to the central NLO prediction 
in the middle panels. The NLO distributions display a 
much smaller dependence on the unphysical renormaliza- 
tion and factorization scales. For our central scale choice, 
the distributions for the first three leading jets soften no- 
ticeably from LO to NLO, while the fourth-jet distribu- 
tion is virtually unchanged. The NLO corrections to the 
behavior of Z, 7* + 4-jet and W + 4-jet production are 
quite similar in this respect [111 ]. 

The bottom panels in fig. [3] show the ratio of Z/W + 
and Z/W~ production both at LO and at NLO. The 
Z/W~ ratio rises with rising pr while the Z/W + ratio is 
roughly flat. Both ratios reflect the rising dominance of 
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the u quark distribution over the d quark with increas- 
ing parton fraction x. Because the Z has an appreciable 
coupling to an initial u quark (unlike the W~ ) , the shape 
of the pt distribution follows more closely the W + case 
than the W~ case, which has a d(x)/u(x) relative sup- 
pression. The excellent agreement between LO and NLO 
ratios for Z/W^ production shows that these ratios are 
under solid perturbative control. 

A comparison of parton-level results to experimental 
data requires estimating the size of non-perturbative ef- 
fects, such as those induced by the underlying event or 
by fragmentation and hadronization of the outgoing par- 
tons. Standard LO parton-shower Monte Carlo programs 
can provide these estimates. As NLO parton-shower pro- 
grams are developed [3l[ , they can use virtual corrections 
computed with BlackHat. We expect non-perturbative 
effects to largely cancel in the Z/W^ ratios. 

In the present study of the Z, 7* + 4-jet process, we 
have imposed cuts typical of Standard-Model measure- 
ments at the LHC. The same code can be used to study 
the size of QCD corrections for observables under cuts 
used in new-physics searches. This will allow the study 
of backgrounds to missing energy signals of new physics, 
arising when a Z boson decays to a pair of neutrinos. 
Ratios such as the Z/W + jets ratios offer highly-reliable 
theoretical predictions. Applying BlackHat along with 
SHERPA brings an unprecedented level of theoretical 
precision to Standard-Model backgrounds, aiding in the 
hunt for new-physics signals at the LHC. 
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